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ABSTRACT 


A general  method  for  representing  the  flow  properties  in  the  three- 
dimensional  boundary  layers  around  ship  hulls  of  arbitrary  shape  is  described. 

It  makes  use  of  an  efficient  two-point  finite-difference  scheme  to  solve  the 
boundary-layer  equations  and  includes  an  algebraic  eddy- viscosity  representa- 
tion of  the  Reynolds-stress  tensor.  The  numerical  method  contains  novel  and 
desirable  features  and  allows  the  calculation  of  flows  in  which  the  circum- 
ferential velocity  component  contains  regions  of  flow  reversal  across  the 
boundary  layer.  The  inviscid  pressure  distribution  is  determined  with  the 
Douglas-Neumann  method  which,  if  necessary,  can  conveniently  allow  for  the 
boundary- layer  displacement  surface.  To  allow  its  application  to  ships,  and 
particularly  to  those  with  double-elliptic  and  flat-bottomed  hulls,  a non- 
orthogonal  coordinate  system  has  been  developed  and  is  shown  to  be  economical, 
precise  and  comparatively  easy  to  use.  Present  calculations  relate  to  zero 
Froude  number  but  they  can  be  extended  to  include  the  effects  of  a water  wave 
and  local  regions  of  flow  separation  which  may  s,tem  from  bulbous-bow  geometries. 
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I.  INTRODUCTION 


A general  method  for  determining  the  local  flow  properties  and  the  over- 
all drag  on  ship  hulls  is  very  desirable  and  particulary  so  with  the  present 
need  to  conserve  energy  resources.  It  is  difficult  to  achieve  for  a number 
of  reasons  including  the  turbulent  nature  of  the  three-dimensional  boundary 
layer,  the  complexity  and  wide  range  of  geometrical  configurations  employed, 
the  possibility  of  local  regions  of  separated  flow  and  the  existence  of  the 
free  surface.  In  addition,  and  although  these  difficulties  may  be  overcome 
in  total  or  in  part,  the  resulting  calculation  method  must  have  the  essential 
features  of  generality,  efficiency  and  accuracy. 

The  purpose  of  this  report  is  to  describe  a general  method  which  is 

capable  of  representing  the  flow  properties  in  the  boundary  layer  around 

ship  hulls  of  arbitrary  shape.  It  is  based  on  the  general  method  of  Cebeci, 

Kaups  and  Ramsey\  developed  for  calculating  three-dimensional,  compressible 

laminar  and  turbulent  boundary  layers  on  arbitrary  wings  and  previously 

proved  to  satisfy  the  requirements  of  numerical  economy  and  precision.  To 

allow  its  application  to  ships  in  general,  and  to  double-elliptic  and  flat- 

bottomed  hulls  in  particular,  an  appropriate  coordinate  system  has  been 

developed.  Previously  described  coordinate  systems,  for  example  a streamline 

2 

system  such  as  that  of  Lin  and  Hall  or  the  orthogonal  arrangement  of  Mi  1 oh 
3 

and  Patel  are  limited  in  their  applicability  and  the  present  nonorthogonal 
arrangement  is  similar  to  that  of  Cebeci,  Kaups  and  Ramsey^. 

The  numerical  procedure  for  solv'nq  the  three-dimensional  boundary-layer 

4 

equations  makes  use  of  Keller's  two-point  finite-difference  method  and  Cebeci 

and  Stewartson's  procedure  in  computing  flows  in  which  the  transverse  velocity 

component  contains  regions  of  reverse  flow.  This  is  in  contrast  to  previous 

2 5 

investigations,  for  example  those  of  Lin  and  Hall  and  Gadd  which  are  limited 
either  to  zero  crossflow  or  to  a unidirectional  and  small  crossflow.  It  is 
also  in  contrast  to  the  previous  methods  of  Chang  and  Patel6  and  Cebeci  and 
Chang'7  which  did  not  have  a good  and  reliable  procedure  for  commuting  the  flow 
in  which  the  transverse  velocity  component  contained  flow  reversal. 


1 


In  representing  turbulent  flow  by  time-averaged  equations,  a turbulence 

model  is  required  and  an  algebraic  eddy-viscosity  formulation,  similar  to  that 

of  Cebeci,  Kaups  and  Ramsey^ , is  used.  This  is  in  contrast  to  the  two- 

£ 

equation  approach  which  Rastogi  and  Rodi  have  applied  to  three-dimensional 

boundary  layers  and  which,  in  principle,  should  be  better  able  to  represent 

flows  which  are  far  from  equilibrium.  The  previous  comparisons  presented  in 
9 10 

Cebeci  ’ demonstrated  that  the  present  eddy-viscosity  model  allows  excellent 
agreement  between  measurements  and  calculations  but  did  not  include  comparison 
with  the  three-dimensional  boundary-layer  measurements  of  Vermeulen^.  Since 
this  data  includes  a strongly  adverse-pressure  gradient  case  which  allows  a 
stringent  test  of  the  present  model,  corresponding  calculations  and  compari- 
sons are  reported. 

The  calculation  method  is  described  in  detail  in  the  following  section 
which  states  the  three-dimensional  boundary-layer  equations  in  curvilinear 
nonorthogonal  coordinates  and  describes  and  discusses  the  required  initial 
conditions,  turbulence  model  and  transformations  in  separate  subsections. 
Section  3 is  devoted  to  the  coordinate  system  which  is  an  essential  feature 
of  the  present  method.  The  numerical  method  is  discussed  briefly  in  section  4 
and  calculated  results  are  presented  in  section  5 which  includes  comparisons 
with  the  measurements  of  Vermeulen^  and  demonstrations  of  the  ability  of  the 
method  to  represent  the  geometry  of  different  hull  configurations  and  to  result 
in  realistic  velocity  and  drag  characteristics.  Summary  conclusions  are  pre- 
sented in  section  6.  Finally,  in  Appendix  A we  present  a description  of  the 
computer  program. 
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II.  BASIC  EQUATIONS 


2.1  The  Boundary-Layer  Equations 

The  governing  boundary-layer  equations  for  three-dimensional  incompres- 
sible laminar  and  turbulent  flows  in  a curvilinear  nonorthogonal  coordinate 
system  are  given  by: 

Continuity  Equation 

~ (uh?  sino<)  + ~ (whj  sino)  + ~ (vhjh2  sino)  = 0 (2.1) 

x-Momentum  Equation 


u au  .w  au  , au 

fi/ax  + fi^az‘+vay" 

_ _ c s c ^ o a_  /P\ 
‘ ~ FTj  ax'  \p  / 

z -Momentum  Equation 
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+ — stf 
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(2.3) 


and  h?  are  the  metric  coefficients  and  they  are,  in  qeneral,  ^unc- 
x and  z;  that  is. 


h-|  = h.j  (x,z) ; h2  = h2(x,z)  (2.4) 

Also,  o represents  the  angle  between  the  coordinates  x and  z.  The  param- 


eters  and 

z = const  and 

k2 

x = 

are  known  as  the  geodesic  curvatures  of  the  curves 
const,  respectively.  They  are  given  by 
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The  parameters  K19  and  K?.  are  defined  by 


M2  sine 


(21  = sine  f \ 2 + HI  az 


reduce  to 


K,  = - 


K12  = ' K1 


1 30  ’ 

^1  3X, 

) + cose  (k2  + TiJ  If)] 

(2.6a) 

1 30' 

) + cose  (K1  + F7I7)] 

(2.6b) 

= ir/2 

and  the  parameters  , Kp,  K^p 

and  K,.| , 

9hl 

3Z 

1 3h2 

• K2  " ‘ 3x 

(2.7) 

K21  K2 

(2.8) 

At  the  edge  of  the  boundary  layer,  (2.2)  and  (2.3)  reduce  to 
u.  3u  w 3u 
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The  boundary  conditions  for  equations  (2.1)  and  (2.3)  are: 


(2.9) 


(2.10) 


y = 0 : u,v,w  = 0 

y = 5 : u = ue(x,z),  w = wg(x,z) 


(2.11a) 

(2.11b) 


2.2  Initial  Conditions  • 

The  solution  of  the  system  given  (2.1)  to  (2.3),  subject  to  (2.11), 
requires  initial  conditions  on  two  planes  intersecting  the  body  along  coordinate 
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lines.  In  general,  the  construction  of  these  initial  conditions  for  three- 
dimensional  flows  on  arbitrary  bodies  such  as  ship  hulls  is  difficult  due  to 
the  variety  of  bow  shapes  which  may  be  extensive  and  complicated.  For  this 
reason,  assumptions  are  necessary  for  them  in  order  to  start  the  calculations. 

In  our  study  we  choose  the  inviscid  dividing  streamline  on  which 
3p/3z  =0  to  be  one  of  the  initial  data  line  (see  figure  1).  In  the  case 


Figure  1.  The  nonorthogonal  coordinate  system  and  the  initial  data  lines 
for  the  ship  hul 1 . 

of  a rectilinear  motion  of  a ship,  this  streamline  runs  along  the  plane  of 
symmetry.  Because  of  symmetry  conditions  w and  ap/8z  are  zero  on  this 
line  causing  (2.3)  to  become  singular.  However,  differentiation  with  respect 
to  z yields  a nonsingular  equation.  After  performing  the  necessary  differ- 
entiation for  the  z-momentum  equation  and  taking  advantaca  of  appropriate 
symmetry  conditions,  we  can  write  the  so-called  longitudinal  attachment-line 
equations  as: 

Continuity  Equation 

fj  (uh2  sine)  + h]  sinew2  + fy  (v  h]h2  sine)  =0  (2.12) 
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x-Momentum  Equation 
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z-Momentum  Equation 


3W 

U 2 

TTj"  3x 


+ V 


aw 
2 

3y 


+ K, 


21 


uw 


3W 


W 


’ K7  -#■ " if  + k21  Vze' 


h{ 


3Wz 


(wVr)z) 


(2.14) 


where  wz  = 3w/3z  and  (wTvr)z  = 3/3z(w'v' ).  These  equations  are  subject 
to  the  following  boundary  conditions: 

y=0:  u=v=wz=0  (2.15a) 

y = 6 : u = t,e’  wz  = wze  (2.15b) 


The  other  initial  data  should  be  selected  near  the  bow  of  the  ship  along 
some  x = const,  coordinate  line  (see  figure  1).  However,  because  of 
the  variety  of  possible  bow  shapes,  approximations  are  necessary.  For  a 
simple,  smooth  bow  section,  where  curvatures  are  small  and  no  separation  is 
expected,  the  flow  along  the  initial  line  can  be  successfully  assumed  to  be 
two-dimensional  without  pressure  gradient,  and  the  governing  two-dimensional 
equations  for  a flat  plate  are  solved.  However,  for  most  general  merchant 
ships,  the  bow  section  is  complicated  and  flow  separation  and  reattachment 
are  expected  because  of  large  curvature  variations  and  adverse  pressure 
gradients;  as  a consequence,  the  boundary-layer  calculations  can  only 
be  performed  downstream  of  the  reattachment  line  (or  point)  where  turbulent 
flow  is  presumed  (since  it  is  unlikely  that  the  flow  remains  laminar  after 
separation  and  reattachment  with  high  Reynolds  number).  Generation  of  the 
initial  data  for  turbulent  flows  is  much  more  involved  if  there  are  no 
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experimental  data  available.  It  requires  sound  mathematical  and  physical 
judgment  and  tedious  trial -and-error  efforts.  We  shall  discuss  this  aspect 
of  the  problem  later  in  the  report. 


In  the  outer  region  e is  defined  by  the  following  formula 


e = 0.0168 
m 


oo 

/ (ute  - ut] 


(2.19) 


where 


2.3  Turbulence  Model 

For  turbulent  flows,  it  is  necessary  to  make  closure  assumptions  for  the 
Reynolds  stresses,  -puTvr  and  -pvTwT.  In  our  study,  we  satisfy  the  require- 
ment by  using  the  eddy-viscosity  concept  and  relate  the  Reynolds  stresses  to 
the  mean  velocity  profiles  by 


-uv  = E, 


m ay  * 


-v  w'  = e. 


(2.16) 


We  use  the  eddy-viscosity  formulation  of  Cebeci  , and  define  e by  two 


separate  formulas.  In  the  inner  region,  em  is  defined  as 
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where 


L = 0.4  y [1  - exp(-y/A)] 
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(2.17) 


(2.18a) 


(2.18b) 


(2.18c) 


ute  = (ue  + we  + 2uewe  cose^ 


(2.20a) 


T 


ut  = ^u2  + w + c°se)  * (2.20b) 

The  inner  and  outer  regions  are  established  by  the  continuity  of  the  eddy- 
viscosity  formula. 

2.4  Transformation  of  the  Basic  Equations 

The  boundary- layer  equations  can  be  solved  either  in  physical  coordinates 
or  in  transformed  coordinates.  Each  coordinate  system  has  its  own  advantage. 
In  three-dimensional  flows,  the  computer  time  and  storage  required  is  an 
important  factor.  The  transformed  coordinates  are  then  favored  because  the 
coordinates  allow  larqer  steps  to  be  taken  in  the  longitudinal  and  transverse 
directions. 

We  define  the  transformed  coordinates  by 


and  introduce  a two-component  vector  potential  such  that 
uh2  sine  = |^-  , wh1  sine  = 
vh1h2  sine  = - (|J  + ff) 

where  \j>  and  41  are  defined  by 

* * (vs]  ue)!^  h2  iinef(x.z.n)  (2.23a) 

* = (vs,  ue)!*  uref/ue  h1  sine  g(x,z,n)  (2.23b) 

and  uref  is  some  reference  velocity. 

Using  these  transformations  and  the  relations  given  by  (2.9),  (2.10)  and 
(2.11),  we  can  write  the  x-momentum  and  z-momentum  equations  for  the  general 
case  as 


(2.22a) 

(2.22b) 
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x-Momentum 

(bf")’  + m^f"  - m2(f')2  - m^f ' g ' + mgf"g  - nig(g')2  + 

af ' 3f \ . /,  3fl  £ ii  33.\ 

" mio V ax  " f ax/  m7  (g  az  ’ f azj 

z-Momentum 

(bg")'  + n^fg"  - m4f'g'  - m3(g')2  + mg  gg"  - mg(f')2  + m12 


(2.24) 


= m 


('•  If  - 5"  %}  ”7(5'  If  - 9"  If) 


10\  ax 

and  their  boundary  conditions  as 

n = 0:  f = f = g - g'  • 0 ; n = n„,  f - 1.  S'  - "e/Vf 

Here  primes  denote  differentiation  with  respect  to  n.  and 

e„ 


u_ 

u. 
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m. 


uref 

ml 

to  r 

S1 

au  \ 

e | 

¥1 

ax  J 

. 1 . + + _ ‘m 

b = I + c , e - — 
m * m v 


The  coefficients  m-j  to  are  given  by 

s 


1 (h,  sine) 

sine  ax  2 


S1  Sue 


m2  " liar  5x~  ■ S1K1  cote  • m3  ■ 'sIK2  u 
1 e 


ref 


cote 


m4  = *1K21  ’ 


_ S1  uref  8ue  . , v uref 

5 " K ~^7  1 K12  u 

2 e e 


m.  = 


6 h^2  sine 


. S1  “ref 

7 " 1I2  ~“7 


u 

\ 


m8  “ s1K2 


. . Uref\ 

<*1  hl  sine  H^J 

W 


csce 


(2.25) 

(2.26) 

(2.27) 


(2.28) 


m 


9 = s-|K-|  csce 
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'10  ~ h 
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ml  1 = m2  + m5  u 


+ mr 


ref 
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» 

v. 


m12  = m4 


w / w \2 

v7  * mfc) 


+ m„  + 


mio  »(.V.  % 

9 ' “ref  8,1  “£f  ,l 


To  transform  the  longitudinal  attachment-line  flow  equations  and  the 
boundary  conditions,  we  use  the  transformed  coordinates  given  by  (2.21)  and 
define  the  two-component  vector  potential  by 


uh2  sine  = §|  , w^  sine  = §* 


vh^h^  sine  = 


(§H 


(2.29) 


with  and  p still  given  by  (2.23).  With  these  variables,  the  longitudinal 
attachment-line  equations  in  the  transformed  coordinates  can  be  written  as 


(bf 


* ”lff"  - m2«”>2  + * ”11  • iff  If  - f"  If)  (2-30) 


(bg")1  + m^'f  - ni^f'g'  - m3(g')2  + m6gg"  - mg(f')2  + 


"12 


= !i  /f  • Ml  _ n» 

h1  \ 3x  9 3x ) 


(2.31) 


The  boundary  conditions  and  the  coefficients  m^  to  m-^  are  the  same  as  in 
(2.26)  and  in  (2.28)  except  now 


n = n_ : 


ze 


ref 


s,  u r 
- 1 ref 


3 ' ^ ue 


m6  = m3 


mg  = 0 , 


'2  3 \uref/  « uref  *7  “ref  3X 


(2.32) 
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ze 


In  terms  of  the  transformed  variables,  the  algebraic  eddy-viscosity 
formulas  as  given  by  (2.17)  to  (2.20)  become 


1 2 


(t”)(  ' /r n 
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1 - exp  I 


(9” 


)%  2 cose  f"g" 


(2.33) 
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III.  COORDINATE  SYSTEM 


Since,  in  general,  a ship  hull  is  a complicated  nondevelopable  surface, 
a Cartesian  coordinate  system  is  not  suitable  for  boundary-layer  calculations. 
Most  existing  merchant  and  naval  vessels  possess  the  following  features:  a 
flat  bottom  [y  = f(x,z)  is  not  a single-valued  function];  a bottom  which 
is  not  parallel  to  the  water  surface;  and  a bow  which  has  a submerged  bulb 
extending  toward  the  origin.  In  addition,  the  problem  is  further  complicated 
by  the  existence  of  a free  surface,  corresponding  to  the  water  level  of  a 
partly-submerged  hull.  The  chosen  coordinate  system  must  be  sufficiently 
general  to  allow  these  various  features  to  be  represented  In  the  boundary- 
layer  calculations. 

The  streamline  coordinate  system  is  superficially  attractive  but  the 
determination  of  the  streamlines,  the  orthogonal  lines,  and  the  associated 
geometrical  parameters  requires  considerable  effort.  They  are  dependent  on 
the  Froude  number,  and  also  on  the  Reynolds  number  if  the  displacement  effect 
is  taken  into  account.  Consequently,  and  in  addition  to  being  hard  to  compute, 
this  coordinate  system  becomes  uneconomical  to  use  when  the  effect  of  the 
Froude  number  and  the  Reynolds  number  are  to  be  systematically  examined. 

A desirable  requirement  of  a coordinate  system  for  the  boundary- layer 

3 

calculations,  is  that  it  be  calculated  only  once.  Miloh  and  Patel  proposed 
an  orthogonal  coordinate  system  which  depends  only  on  the  body  geometry  and 
is  calculated  once  and  for  all.  This  coordinate  system  has  been  applied  by 
Chang  and  Patel^  to  boundary-layer  calculations  on  two  simple  ship  hulls: 
ellipsoid  and  double  elliptic  ship.  One  of  the  coordinates  is  taken  as  lines 
of  x = x = constant  and  the  other  as  z(x,z)  = constant,  which  is  ortho- 
gonal to  x = constant  lines  everywhere  on  the  ship  hull,  and  is  obtained 
from  the  solution  of  the  differential  equation 


dz 

dx 


fxfi 


(3.1) 


T+rf 

z 

Here  y = f(x,z)  defines  the  ship  hull,  and  (x,y,z)  denote  the  Cartesian 
coordinates.  The  major  advantage  of  this  coordinate  system  is  its  simplicity. 
Because  one  of  the  coordinates  is  subject  to  the  condition  (3.1),  there  is  no 
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I* 

f. 


guarantee  that  the  boundaries  of  the  ship  hull  are  coincident  with  the  coord- 
inate lines.  Furthermore,  for  a ship  with  flat  bottom  for  which  y = f(x,z) 
is  not  a single-valued  function,  one  of  the  coordinates  cannot  be  calculated 
from  (3.1).  The  coordinate  system  is  limited,  therefore,  to  some  special 
geometries  only. 

In  this  study  we  adopt  a nonorthogonal  coordinate  system  similar  to  that 
developed  by  Cebeci , Kaups  and  Ramsey^  for  arbitrary  wings.  It  is  based  on 
body  geometries  only  and,  hence,  it  is  calculated  once  and  for  all.  In 
addition,  the  system  can  deal  with  the  peculiar  features  of  most  merchant  and 
naval  vessels  discussed  previously.  The  details  of  this  coordinate  system  are 
ucs'^ibed  briefly  in  the  following  paragraph. 

Now  consider  the  ship  hull  as  given  in  the  usual  Cartesian  coordinate 
system;  that  is,  x along  the  ship  axis,  y and  z in  the  crossplane  (see 
figure  1).  We  select  x = x = constant  as  one  of  the  coordinates  and  the 
other  coordinate  z lies  in  the  yz-plane.  Because  the  coordinate  system  is 
nonorthogonal,  we  are  free  to  select  the  values  of  z in  the  plane  to  satisfy 
the  condition  that  the  boundary  lines  of  the  ship  hull  are  coincident  with 
z = constant  coordinate  lines.  There  are  several  ways  of  defining  the  z-values. 
Here  z is  determined  by  mapping  each  yz  crossplane  into  a half  or  full 
unit  circle  depending  on  whether  the  crossplane  intercepts  the  free  surface 
or  is  completely  submerged.  The  polar  angle,  normalized  by  ir  or  2ir  on  the 
unit  circle,  is  taken  as  z-values.  The  z-values  then  range  from  0 to  1 on 
each  crossplane.  The  advantage  of  the  mapping  method  is  that  equi -Interval , 
z=constant  coordinate  lines  are  automatically  concentrated  in  the  region  of 
large  curvature  where  the  boundary-layer  characteristics  are  expected  to 
vary  greatly.  Hence  the  number  of  z=constant  coordinate  lines  can  be  reduced 
without  loss  of  accuracy. 

There  are  several  methods  available  for  the  mapping  of  an  arbitrary  body 

onto  a unit  circle.  Here  we  use  the  numerical  mapping  method  developed  by 
12 

Halsey  . It  makes  full  use  of  Fast  Fourier  Transform  techniques  and  has  no 
restrictions  on  the  shape  of  the  body  to  be  mapped.  To  map  a smooth  cross- 
plane onto  a unit  circle,  the  procedure  is  fairly  easy.  If  there  are  inner 
corner  points,  or  trailing-edge  and  leading-edge  corner  points  (see  figure  2) 
caused  by  the  reflection  of  the  crossplane,  they  must  be  removed  before  mapping 
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Fiqure  2.  Notation  of  corner  points  used  in  the  mappinq  procedure. 

is  performed  to  improve  numerical  accuracy  and  to  provide  rapid  convergence. 

The  inner  corner  points  are  rounded  off  by  using  Fourier  series  expansion 

technique  and  the  leading-edge  and/or  trai ling-edge  corner  points  and  removed 

12 

by  using  the  Karman-Trefftz  mapping.  For  details  see  Halsey  . 

To  use  the  mapping  method  to  find  the  coordinate  system,  it  is  only 
necessary  to  define  the  ship  hull  as  a family  of  points  in  the  x=constant 
planes,  to  locate  the  intersection  of  the  ship  hull  and  the  free  surface, 
and  to  indicate  whether  corner  points  exist.  The  data  in  each  plane  is  then 
mapped  into  a unit  circle  as  y vs  z and  z vs  z and  interpolated  for 
constant  values  of  z.  Another  set  of  spline  fits  in  the  planes  z=constant 
for  y vs  x and  z vs  x completes  the  definition  of  the  coordinate  system. 
The  lines  formed  by  the  intersection  of  the  planes  x=constant  and  z=constant 
with  the  hull  constitute  the  nonorthogonal  coordinate  net  on  the  surface,  and 
the  third  boundary-layer  coordinate  is  taken  as  the  distance  normal  to  the 
surface  in  accordance  with  first-order  boundary-layer  approximation. 

Since  the  spline-fitting  also  yields  derivatives,  the  metric  coefficient 
and  the  geodesic  curvatures  of  the  coordinate  lines  can  be  calculated  from 
the  formulas  given  below. 
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The  metric  coefficients: 


* (£): 

t ■ at  * i( 


The  angle  between  the  coordinate  lines: 


cose  = Kfc  (ll)x(S)z  + (If )z  (If), 


The  geodesic  curvature  of  the  z=constant  line 


/li\  /li\  a2* 

)\3Z  1 \ 3X  / * \3x  / \3z  / \3x  3>r  ’ 3X  3X7 

/L'  ;x'  'z  ' 'z ' 'xJl' 


(3.2a) 


The  geodesic  curvature  of  the  x=constant  line 


K2  = 

2 h^h^  si 


(21\  /32\  /a£\  /3z\  ai  32z\  ,,  0 

ino  \d2L\*x)7  \3X/  13 zMaz  3z2  3Z  3Z7/  VJ,0J 


The  other  parameters  K12  and  K?1  are  calculated  from  (2.6).  It  may  be 
noted  that  K-j  and  1^  can  also  be  obtained  from  (2.5).  This  provides  a 
check  on  the  expressions  given  by  (3.4)  and  (3.5). 

In  boundary-layer  calculations,  we  need  the  inviscid  velocity  compon- 
ents along  the  surface  coordinates.  Let  V be  the  total  velocity  vector  on 
the  hull,  (u,v,w)  the  corresponding  velocity  components  in  the  Cartesian 
coordinates,  and  (ue,  we)  in  the  adopted  surface  coordinates.  As  can  be 
seen  from  figure  3, 


V • t-|  — cose  V • t2 
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(3.7) 
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Here  and  1 2 are  the  unit  tangent  vectors  along  x ar.d  z coordinates 
and  are  given  by 


*7 


i + 


(3.8) 


(3.9) 


-V 

With  the  definition  of  V and  with  the  use  of  (3.8)  and  (3.9),  equations 
(3.6)  and  (3.7)  can  be  written  as 


u 


e 


- (sz  \ 


1 cose 
Ti? 


(3.10) 


IV.  NUMERICAL  METHOD 

We  use  the  Box  method  to  solve  the  boundary-layer  equations  given  in 

Section  II.  This  is  a two-point  finite-difference  method  developed  by  Keller 

and  Cebeci . This  method  has  been  applied  to  two-dimensional  flows  as  well  as 

three-dimensional  flows  and  has  been  found  to  be  efficient  and  accurate. 

Descriptions  of  this  method  have  been  presented  in  a series  of  papers  and 

reports  and  a detailed  presentation  is  contained  in  a recent  book  by  Cebeci 
1 3 

and  Bradshaw  . 

In  using  this  numerical  method,  or  any  other  method,  care  must  be  taken 
in  obtaining  solutions  of  the  equations  when  the  transverse  velocity  compon- 
ent w contains  regions  of  flow  reversal.  Such  change  in  w-profiles  will 
lead  to  numerical  instabilities  resulting  from  integration  opposed  to  the 
flow  direction  unless  appropriate  changes  are  made  in  the  integration  pro- 
cedure. Here  we  use  the  procedure  developed  by  Cebeci  and  Stewartson.  In 
this  new  and  very  powerful  procedure,  which  follows  the  characteristics  of 
the  locally  plane  flow,  the  direction  of  w at  each  grid  point  across  the 
boundary  layer  is  checked  and  difference  equations  are  written  accordingly. 

At  each  point  to  be  calculated,  the  backward  characteristics  which  determine 
the  domain  of  dependence,  are  computed  from  the  local  values  of  the  velocity. 
Since  the  characteristic  must  be  determined  as  part  of  the  solution,  in  the 
calculation  procedure  a Newton  iteration  process  is  used  to  correctly  deter- 
mine the  exact  shape  of  the  domain  of  dependence. 

To  illustrate  the  basic  numerical  method,  we  shall  at  first  consider  the 
solution  of  the  longitudinal  attachment-line  equations  (2.30)  and  (2.31)  and 
then  the  solution  of  the  full  three-dimensional  flow  equations  are  given  by 
(2.24)  and  (2.25).  We  shall  not  discuss  the  Cebeci -Stewartson  procedure  for 
computing  three-dimensional  flows  with  the  transverse  velocity,  w,  con- 
taining flow  reversal  since  that  procedure  will  be  fully  described  in  a forth- 
coming paper. 


4 . 1 Difference  Equations  for  the  Longitudinal  Attachment-Line  Equations 

According  to  the  Box  method,  we  first  reduce  the  equations  (2.30),  (2.31), 
(2.32)  and  (2.26)  into  a system  of  five  first-order  equations  by  introducing 
new  dependent  variables  u(x,z,n),  v(x,z,n),  w(x,z,n),  t(x,z,n)  and  e(x,z,n). 
Equations  (2.30)  and  (2.31)  can  then  be  written  as 


u'  = V 


w‘  = t 

2 


(bv) 1 + uv  - + 

2 


e = m^u  + nigW  + m 

The  boundary  conditions  (26)  and  (32)  become 

n=0:  u = w = 0 = 0 

n = n„,  : u = 1 , w = wze/uref 


(4.1a) 

(4.1b) 

3U 

"’ll  " mi0  U 3X 

(4.1c) 

2 9W 

-m9u  + *12  “ m10  U SI 

(4. Id) 

i + min  — 

10  ax 

(4 . 1 e ) 

(4.2) 


We  next  consider  the  net  rectangle  shown  in  figure  4 and  denote  the  net 
points  by 


xo-0 
"o  ” 0 


x = x . + k 
n n-1  n 


nj  = n 


. , + h . 
j-1  J 


n = 1 , 2,  N 

j = 1.  2,  J 


ni-i/2 

Vi 


(*  ^n— 1/2 
(» 


Figure  4.  Net  rectangle  for  the  longitudinal  attachment-line  equations. 

We  approximate  the  quantities  (u,v,g,t,o)  at  points  (xn,nj)  of  the 
net  by  funtions  denoted  by  (uj,  vj * wj»  • ej) • We  also  emPloy  the  notation 
for  points  and  quantities  midway  between  net  points  and  for  any.  net  function  s" 
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(4.3) 


Vi/2  S?K*  Vi] 


'Tw-z  "j"!'1' 


"j-l/2  = 2 (nj  + "j-11 
Sj-l/2  ' 2 (sj  sj-l* 


The  difference  equations  which  are  to  approximate  (4.1)  are  formulated  by 
considering  one  mesh  rectangle  as  in  figure  4.  We  approximate  (4 . 1 a ,b ) 
using  centered  difference  quotients  and  average  them  about  the  midpoint 

( xn » n j _i /2 ^ of  the  PlP2‘ 


,-l,  n n v n 
hj  (uj  Uj-1)  ' j-l/2 

hj  (wj  Wj-1}  ' j-l/2 


(4.4a) 

(4.4b) 


Similarly,  (4.1c,d,e)  are  approximated  by  centering  them  about  the  midpoint 
( xn - 1 / 2 * n i - 1 / 2 ) the  rectan9le  P1 P2P3P4 * This  g-ives 

hj'tih.);  - (fc.jj.,1  * <«>j_1/2  - (m2  ♦ °nH»%y2  - -*?,  «.w 


hj  * (0t>j-l/2  (m4  + °n>(uw>"-l/2  m3*"Z)j-l/2 

- m9 <u^> j-l/2  * “r.[“j-V2Uj-l/2  _ "j-l/2Wj-l/2^  ' Sj-\/2  ~ ”l2 


"S’Cj  - 9J-1  > - K + 2“n>b". 


- mnwn  - Tn*l 
1/2  m6W j -1/2  " j-l/2 


(4.4d) 
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Here 
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Kj’Khv)-'  - (bv)J;J]  ♦ (ev)-J/2 


n-1  ,..2xn-l 


+ m 
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m2  ^ j-l/2  ' "'ll 


(4.5a) 
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vn-1 
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Tn-1  _ , n-1 

Jj-l/2  ~ ~2nnuj -1/2 

m n'1/2 
10 

xn  “ xn-l 


°j=!  > 


n-l  n-1  n-1  n-1  | 

m1  uj-l/2  m6  wj-l/2) 


(4.5c) 


(4.5d) 


4.2  Difference  Equations  for  the  Full  Three-Dimensional  Equations 

The  difference  equations  for  the  full  three-dimensional  equations,  as 
given  by  (2.24)  and  (2.25),  are  again  expressed  in  terms  of  a first-order 
system.  With  the  definitions  given  by  (4.1a)  and  (4.1b),  they  are  written 


2 2 3 u 3 u 

(bv)1  + ov  — nipU  — m^uw  — iiigiv  + u — + m^  w ~ 

22  3w  )w 

(b t ) * + ot  — m^uw  — — ntgii  + in^  = u * m7  w 

o'  = m.u  + m,w  + “ + in-,  ^ 

1 6 10  dx  7 ;iz 

Their  boundary  conditions,  (2.26)  become: 


n = 0: 

n - nj  u = 1 , 


u = w = o =0 

W ■ VVf 


(4.6a) 


(4.6b) 


(4.6c) 


(4.7a) 

(4.7b) 


The  difference  equations  for  (4.1a)  and  (4.1b)  are  the  same  as  those 
given  by  (4.4a)  and  (4.4b):  they  are  written  for  the  midpoint  (xn, 
n j_i /2 ) the  net  cube  shown  in  figure  5;  that  is 

h-1fun.i  = vn»i  h”^(wn’^  — wn,i)  = tn’^  (4  8) 

j (Uj  j-r  Vj-l/2 ’ hj  (wj  wj-l'  j-1/2  '''‘''I 

The  difference  equations  whcih  are  to  approximate  (4.6a,b,c)  are  rather 
lengthy.  To  illustrate  the  difference  equations  for  these  three  equations, 
we  consider  the  following  model  equation 


(bv)'  ♦ ev  v u |£  + «,»  f»- 


(4.9) 


J 


'JW  . 


/ 


Figure  5. 


Net  cube  for  the  difference  equations  for  three-dimensional 
flows,  Wj  > 0. 


The  difference  equations  for  (4.9)  are: 

hj1[(E7)j  “ (57)j-l]  + ^j-1/2  + (ml  1 ) j-l/l 

= (m10)i-l/2  uj-l/2  + m7  i-1/2  W j -1/2  r. 
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Here,  for  example 
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“n  = T (“j*1  + “j’1'1  + u";5  + “"-l"11 

S'1,1  tu":i 
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z_  = 0 z;  = z.  , + r,  i = 1 , 2,  . . . , I 
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(4. 1C) 


(4.11) 


(4.12) 


4.3  Solution  of  the  Difference  Equations 

The  difference  equations  (4.4)  for  the  longitudinal  attachment-line  flow 
and  the  difference  equations  for  (4.6)  are  nonlinear  algebraic  equations.  We 
use  Newton's  method  to  linearize  them  and  then  solve  the  resulting  linear 
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system  by  the  block-elimination  method  discussed  by  Keller  . A brief  descrip- 
tion of  it  will  be  given  for  the  longitudinal  attachment-line  equations. 

Using  Newton's  method,  the  linearized  difference  equations  for  the  system 
given  by  (4.4)  are: 

6Uj  — jjJ-Uvj  + -j)  = (r^j  (4.13a) 

h . 

6Wj  - 6Wj_.j  -^'('Stj  + 4tj_-j)  = (r2)j  (4.13b) 

<'l>jsvJ  * (42>s,j-l  * * (t4>jS9M  + <t5)jiuj 

+ (c6)j6Uj _1  = (r3)j  (4.13c) 


♦ UjJjJtj.,  * * (85)jSWj 

+ KV^j-i  + (B7)j6Uj  + ^8VUJ-l  = (r4^j  (4J3d) 

(ct1  )j60j  + (o2)j60j-l  + (a3)jiuj  + (a4)j<Suj-l  + (o5)j6wj 

+ (og)j6Wj_i  = (r5)j  (4.13e) 

Here  we  have  dropped  the  superscripts  n,  i and  have  defined  (rk)^,  (ck)j. 
(6k)j  and  (ok).  by 


rl  = Uj-1  ‘Uj  + hjVj-l/2 


(4.14a) 


r2  = Wj-1  wj  + hjtj-l/2 


(4.14b) 


(r3*j  = Rj-l/2  mll 


C(bv)j-l/2  + (0v>j-l/2  - <m2  + an)(u  >j-l/2]  (4*14c) 


(rVj  ' Sj-l/2  m12  C (bt) j-1/2  + ^ot^j-l/2  K + an^uw)j-l/2 

-m3(w2)j-i/2  ~ m9(u^j-i/2  + an(wj-l/2uj-l/2  ~ uj-l/2wj-l/2)^4J4d) 


= Tj-l/2  [0j-l/2  (ml  + 2an)uj-l/2  m6w  j - 1 / 2 ^ 

(4.1 4e ) 

Ul}j  = hf  + F°j 

J 

(4.15a) 

(4.15b) 

(c3>j  = \ vj 

(4.15c) 

(c4Jj  = \ Vj-1 

(4.1 5d) 

U5}j  = - {m2  + an)uj 

(4.15e) 

(c6)j  = - (m2  + «n)Uj.-| 

(4.1 5f ) 

= 

(4.16a) 

(e2)j  = (c2)j 

(4.16b) 

(e3*j  = \ 

(4.16c) 

= \ lj-l 

(4 . 1 6d ) 

(e5)j  = - ? ^m4  + °VUj  ~~  m3wj  \ anUj-l/2 

(4.1 6e ) 

1 1 n 1 

= ~1  H + an)uj-l  _m3wj-l  ~2  anuj-l/2 

(4 . 1 6f ) 

(B7)j  3 ~ \ (m4  + an)Wj  - m9uj  + \ anwj-l/2 

(4.1 6g ) 

(e8*j  3 ~ \ K + arJwj-l  _m9uj-l  + 2anwj-l/2 

(4.1 6h ) 

(ai}j  a h 

(4.17a) 
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* 2an> 

(4.1 7d) 

(°5>j 

1 

' 2 

m6 

(4. 17e) 

K>j 

1 

2 

m6 

(4.17f ) 

The  boundary  conditions  become 

6UQ  = 6Wq  = <50q  = 0,  6Uj  = <5Wj  = 0 (4.18) 

The  solution  of  the  linear  system  given  by  (4.13)  and  (4.18)  is  obtained 
by  using  the  block  elimination  method.  According  to  this  method,  the  system 
is  written  as 

r. 

Here 
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The  A.,  B.,  C.  and  A denote  5x5  matrices.  The  solution  of  (4.19)  is 
J J J 13 

obtained  by  the  procedure  described  In  Cebeci  and  Bradshaw  . 
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V.  RESULTS 


5.1  Turbulent  Flow  Calculations  for  a Curved  Duct  and  Comparison  with 
Experiment 

The  turbulence  model  described  in  Section  II  has  been  used  with  consider- 
able success  to  compute  a wide  range  of  two-dimensional  turbulent  boundary 
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layers,  see  for  example  Cebeci  and  Smith  . The  model  has  also  been  used  to 
compute  three-dimensional  flows  and  again  is  found  to  yield  accurate  results 
(see  for  example  Cebeci^’^  and  Cebeci,  Kaups  and  Moser*^).  To  further  test 
the  model  for  three-dimensional  flows,  we  have  considered  the  experimental 
data  taken  in  a 60°  curved  duct  of  rectangular  cross  section.  Fiqure  6 shows 
a sketch  of  the  flow  geometry.  The  experimental  data  are  due  to  Vermeulen1^. 
Here  z denotes  the  distance  from  the  outer  wall,  measured  along  normals  to 
the  wall;  x denotes  the  arc  length  along  the  outer  wall;  and  y denotes 
distance  normal  to  the  plane  x,z. 


Figure  6.  Coordinate  system  and  notation  for  the  curved  duct. 

To  compare  the  computed  results  with  the  data,  it  is  necessary  to  specify 
the  initial  profiles  given  by  experiment.  This  can  be  done  in  a number  of 
ways.  In  the  study  reported  by  Cebeci,  Kaups  and  Moser^  the  profiles  were 
generated  by  using  Coles'  velocity  profile  formula.  That  formula,  which 
represents  the  experimental  data  rather  well  for  two-dimensional  flows,  was 
not  very  satisfactory  for  three-dimensional  flows.  Here  we  abandon  the  use 
of  Coles'  formula  in  favor  of  Thompson's  two-parameter  velocity  profiles  as 
described  and  improved  by  Galbraith  and  Head^.  According  to  this  formula. 
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the  dimensionless  u/ug  velocity  profile  is  given  by 


— = y / “ — \ ♦ 0 -YJ 

ue  sVueL_ 


(5.1) 


Here  y is  an  intermittency  factor  defined  by  the  following  empirical 
formulas: 


0 < < 0.05 

5o  " 


0.05  < f-  •<  0.3 
0 — 

0 


0.3  < f-  < 0.7 
6o  “ 


0.7  < < 0.95 

6o  ” 


YS  - 1 


ys  = 1 - 2.64214  - 0.05j 

= 4.4053^  - 0.5  j - 1 . 8502 ^ 0.5^+  0.5 

/ \2 
= 2.64214  - 0.05  j 


(5.2) 


y = 2.64214 
's 


X-  > 0.95 
6o 


Ys  = 0.0 


The  dimensionless  velocity  profile  for  the  inner  layer,  that  is,  (u/ue\-nner’ 
is  given  by 


+ . 
y < 4 


+ + 
u = y 


4 < y+  < 30  u+  = c1  + c2  In  y+  + c3(ln  y+)2  + c4(ln  y+)3 

y+  > 30  u+  = 5.50  In  y +5.45 


(5.3) 


Here  c,  = 4.187,  c9  = -5.745,  c.,  = 5.110,  c.  = -0.767,  y = yu  /v, 

I «i  jf\  + C 3 H 

u = (t  /p)  , u = u/u  and  is  a parameter  which  is  a function  of 

t w t o 

9,  c ^ and  H. 

To  find  the  functional  relationship  between  6Q,  c^,  e and  H,  we  use 
the  definitions  of  displacement  thickness  6*  and  momentum  thickness  e. 
Substituting  (5.1)  into  the  definition  of  6*.  after  some  algebra,  we  get 


6^  (]  “ R^*)  = °-6+^P  A4  ln  A3  *2  - M7) 


(5.4) 


where 

A1  = 50.679,  A2  = 1.1942,  A3  = 0.7943,  A4  = 1.195 

An  expression  similar  to  that  given  by  (5.4)  can  also  be  obtained  if  we 
substitute  (5.1)  into  the  definition  of  e.  However,  the  resulting  expression 
is  quite  complicated.  For  this  reason,  the  expression  for  e is  obtained 
numerically,  and  for  a given  value  of  e and  H,  the  corresponding  values 
of  c^  and  6Q  are  computed  from  that  equation  and  from  (5.4). 

Equation  (5.1)  is  recommended  for  two-dimensional  flows.  Here  we  assume 
that  it  also  applies  to  the  streamwise  velocity  profile  by  replacing  u/ug  by 
us/uSe  with  cf  now  representing  the  streamwise  skin-friction  coefficient. 

In  order  to  generate  the  crossflow  velocity  component  (un/uSe),  we 
use  Mager's  expression  and  define  un/uSg  by 


with  the  limiting  crossflow  angle  8^  obtained  from  the  experimental  data. 

Once  the  streamwise  and  crossflow  velocity  profiles  are  calculated  by 
the  above  procedure,  we  compute  the  velocity  profiles  u/ug  and  w/wg  in 
the  orthogonal  directions  x and  z by  the  following  relationships 


u^ 

u 


Jse 


u„  w 
_n 

uc  u 
Se  ue 


w_ 

w„ 


s_+  VA 
Se  use  we 


(5.6a) 

(5.6b) 


Figure  7 shows  a comparison  of  generated  and  experimental  total  velocity 
profiles  along  the  line  A.  As  can  be  seen,  the  above  discussed  procedure 
for  generating  the  initial  velocity  profiles  from  the  experimental  data  is 
quite  good.  This  is  important  for  an  accurate  evaluation  of  a turbulent 
model  especially  for  three-dimensional  flows.  Here 
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Figure  7.  Comparison  of  generated  initial  total  velocity  profiles  with 
Vermeulen's  data. 


(5.7 


The  solution  of  the  boundary-layer  equations  also  requires  the  specifica 
tion  of  the  metric  coefficients  and  the  geodesic  curvatures.  They  are  calcu- 


lated from  the  following  expression: 


h2  = 1.0, 


straight  section 
curved  section 

K2  = 0 (5.8 


K 


1 


0 


straight  section 


1/(R  — z)  curved  section 


immj 


A comparison  of  calculated  and  experimental  values  of  streamwise  momentum 
thickness  e^,  shape  factor  H^,  skin-friction  coefficient  cf  and  limit- 
ing crossflow  angle  6 is  shown  in  figures  8,  9,  10  and  11,  respectively, 

w 

along  the  lines  B,  C,  D,  E.  Here  the  limiting  crossflow  angle  is  computed 
from 


tane 

W 


w /u  [ (u  Qf/u  )g"  — f"] 
e eLV  ref  e 3w  w 

(w„/u  (u  f/w  )g"  + f" 
' e'  e ref  e aw  w 


(5.9) 


Figure  8.  Comparison  of  computed  Figure  9.  Comparison  of  computed 
momentum  thickness  with  shape  factor  with 

Vermeulen 1 s data.  Vermeulen's  data. 


Figure  10.  Comparison  of  computed  Figure  11.  Comparison  of  computed 
skin  friction  coefficient  limiting  crossflow  angle 

with  Vermeulen's  data.  with  Vermeulen's  data. 
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Figures  12  and  13  show  a comparison  of  calculated  and  experimental  total 
velocity  profiles  and  crossflow  angle  profiles  along  the  lines  C and  E. 
Here  the  crossflow  angle  is  computed  from 


we/ue  [(uref/we)g‘  "f']  (5.10) 

sine  = p 

Q/Qe  (Qe/“er 

As  in  figures  8 through  11,  again  the  agreement  between  calculated  results 
and  experiment  is  very  good.  The  computed  results  follow  the  trend  in  the 
experimental  data  well  and  indicate  that  the  present  turbulence  model,  as  in 
two-dimensional  flows,  is  quite  satisfactory  for  three-dimensional  flows. 


5.2  Results  for  a Double  Elliptic  Ship  Model 

To  test  our  method  for  ship  hulls,  we  have  considered  two  separate  hulls. 
The  first  one,  which  is  discussed  in  this  section,  is  a double  elliptic  ship 
whose  hull  is  given  analytically.  The  second  one,  which  is  discussed  in 
section  C,  is  ship  model  5350  which  has  a rather  complex  shape.  Its  hull  is 
represented  section-by-section  in  tabular  form  and  contains  all  the  features 
of  most  merchant  and  naval  vessels.  It  proves  an  excellent  test  case  to  study 
the  computational  difficulties  associated  with  real  ship  hulls. 

The  double  elliptic  ship  model  can  be  analytically  represented  by 


C 1 

rv> 

1/2 

/;\2 

y = f(x,z)  = B 

i 

n> 

1 -(§)] 

It  has  round  edges  except  for  the  sharp  corners  at  x = +L  and  z = +H.  The 
body  of  L:H:B  = 1.0:0.125:0.1  together  with  the  nonorthogonal  coordinate 
nets  on  the  hull  is  shown  in  figure  14. 

The  potential -flow  solutions  were  obtained  from  the  Douglas-Neumann 
computer  program  for  three-dimensional  flows.  To  get  the  solutions,  120  con- 
trol elements  on  the  surface  were  used,  12  along  the  x-direction  and  10  along 
the  z-di recti  on. 

Before  we  describe  our  boundary-layer  calculations,  it  is  useful  to 
discuss  the  pressure  distribution  for  this  body  shown  in  figure  15.  As  can  be 
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Figure  12.  Comparison  of  computed  total  velocity  profiles  with 
Vermeulen’s  data. 


Figure  13.  Comparison  of  computed  crossflow  angle  with  Vermeulen's 
data. 


Fiqure  14.  Three-dimensional  picture  of  double  elliptic  ship  model  with 
the  nonorthogonal  coordinate  system. 
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seen  from  the  figure*  the  longitudinal 
pressure  gradient  is  initially  favorable 
in  the  bow  region  and  then  adverse  up 
to  the  midpoint  of  the  body.  This  is 
followed  by  a region  of  favorable 
pressure  gradient  and  then  by  a sharp 
adverse  pressure  gradient  very  close  to 
the  stern.  The  transverse  pressure 
gradient  varies  in  a more  complex 
manner.  Near  the  bow  the  pressure 
decreases  from  the  water  surface  to 
9 minimum  and  then  increases  as  the 
keel  is  reached.  As  the  flow  moves 
downstream,  the  location  of  the  mini- 
mum pressure  moves  up  and  reaches  the 
water  surface  at  about  x/L  = -0.80. 

The  minimum  pressure  remains  at  the 
water  surface  to  about  x/L  = 0.80  and  then  moves  toward  the  keel.  As  a 
result,  near  the  bow  and  the  stern,  one  may  expect  flow  reversal  of  the 
transverse  profiles  while  on  the  remaining  major  part  of  the  ship  hull,  the 
transverse  across  the  boundary  layer  does  not  reverse  direction  from  the  keel 
to  the  water  surface.  This  conclusion  is  drawn  from  considering  the  pressure 
gradients  only.  The  real  situation  may  be  somewhat  modified  because,  in 
addition,  there  are  the  upstream  effects  and  the  curvature  effects  on  the 
flow  characteristics. 

The  boundary-layer  computation  starts  with  turbulent  flow  from  x/L  = -0.90. 
We  have  tried  to  start  the  computation  from  x/L  = -0.97  and  x/L  = -0.95. 
However,  flow  separation  was  observed  at  x/L  = -0.90  near  the  keel  due  to 
the  sharp  curvature  and  adverse  pressure  gradient  in  the  bow  region  and  can  be 
seen  from  figure  15.  In  the  previous  calculations  of  Chang  and  Patel6  and 
Cebeci  and  Chang^,  the  flow  separation  near  the  bow  was  not  found  due  to  the 
orthogonal  coordinate  system  they  adopted  in  which  the  second  net  point  from 
the  keel  is  so  far  from  the  keel  that  the  region  of  adverse  pressure  gradient 
is  omitted. 


the  double-elliptic  ship. 
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In  our  boundary-layer  calculations,  we  have  used  40  points  along  the 

x-di recti  on  and  16  points  along  the  z-di  recti  on.  In  the  normal  direction, 

we  have  taken  approximately  40  points.  The  nonuniform  grid  structure  des- 
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cribed  in  Cebeci  and  Bradshaw  is  employed  in  the  normal  direction  so  that 
the  grid  points  are  concentrated  near  the  wall  where  the  velocity  gradients 
are  large. 

Some  of  the  computed  results  for  RL  = 107  are  shown  in  figures  16  to 
18.  Figure  16  shows  the  transverse  distributions  of  the  pressure  coefficients, 
cp,  local  skin-friction  coefficient,  cf,  the  shape  factor  H^,  the 
Reynolds  number  based  on  the  momentum  thickness,  R and  the  limiting  cross- 
flow  angle  for  x/L  = -0.85,  -0.50,  0.0,  0.25,  0.50,  0.75.  As  can  be  seen  from 
these  figures,  the  boundary-layer  parameters  vary  greatly  near  the  keel  where 
the  curvatures  and  the  pressure  gradients  are  large  and  remain  almost  unchanged 
near  the  water  surface  where  the  curvatures  and  the  pressure  gradients  are 
small.  Except  at  x/L  = -0.85,  the  limiting  crossflow  angle  is  positive.  This 
implies  that  the  crossflow  near  the  wall  moves  from  the  keel  to  the  free 
surface  as  predicted  from  the  pressure  distribution.  Figure  17  shows  typical 
longitudinal  and  transverse  velocity  profiles  at  z = 0.6  for  several  values 


Figure  16.  Computed  cp,  cf,  Hn,  Re  and  &.  for  the  double-elliptic  ship 
model  for  R.  = 107at  (a)  x/L  = -0.85,  (b)  x/L  = -0.50,  (c) 

x/L  = 0.0,  (d)  x/L  = 0.25,  (e)  x/L  = 0.50,  (f)  x/L  = 0.75. 
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of  x/L  and  figure  18  shows  typical 
transverse  velocity  profiles  at  x/L 
= -0.2  for  several  values  of  z.  As 
can  be  seen  from  figures  17(b)  and  18, 
the  transverse  velocity  component  under- 
goes drastic  changes  in  the  longitudinal 
and  transverse  directions  under  the 
influence  of  pressure  gradient  and  body 
geometry.  As  was  discussed  before, 
when  the  transverse  velocity  changes 
sign  across  the  boundary  layer  and  con- 
tains regions  of  reverse  flow,  numerical 
instabilities  result  from  integration 
opposed  to  flow  direction  unless  appro- 
priate changes  are  made  in  the  integra- 
tion procedure.  The  new  numerical  pro- 
cedure of  Cebeci  and  Stewartson  handles 


Figure  18.  Computed  transverse  veloc- 
ity profiles  for  the  double- 
elliptic  ship  model  at 


x/L  = -0.2. 

this  situation  very  well  and  does  not  show  any  signs  of  breakdown  resulting 
from  flow  reversal  of  transverse  velocity  component. 


5.3  Results  of  Ship  Model  5350 

The  ship  model  5350,  unlike  the  one  discussed  in  section  B,  is  a realistic 
tanker  model.  The  geometry  of  the  hull  is  so  complicated  that  it  is  repre- 
sented in  tabular  form  section  by  section.  The  model  possesses  all  the  special 
features  of  existing  merchant  and  naval  vessels,  that  is,  a bottom  which  is 
flat  and  not  parallel  to  the  still-water  surface  and  an  extended  bow  completely 
submerged  under  the  water  surface,  and  consequently  serves  as  an  excellent  case 
on  which  to  apply  our  method. 

Figure  19  shows  a three-dimensional  picture  of  this  ship  model  together 
with  our  nonorthogonal  coordinate  system.  We  see  from  this  figure  that,  as 
a by-product  of  the  mapping  method  discussed  earlier,  the  z=const.  coordinate 
lines  are  concentrated  in  the  bow  and  corner  regions  where  the  curvature  is 
large.  Figure  20  shows  different  cross-sections  (indicated  by  solid  lines) 
and  interpolated  values  obtained  by  a cubic-spline  method  (indicated  by 
circles)  from  which  the  geometric  parameters  are  obtained. 
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Figure  19.  Three-dimensional  view  of  ship  model  5350  with  the  nonorthogonal 
coordinate  system. 


INPUT  SHIP  FORM 

INTERPOLATION  BY  CUBIC  SPl  INE  FUNCTION 


Figure  20.  Body  plan  for  ship  model  5350. 


The  inviscid  velocity  distribution  for  the  model  is  obtained  by  using  the 
Douglas-Neumann  method  treating  the  model  as  a double  ship  model.  Figure  21 
shows  the  pressure  distribution  for  the  entire  ship  and  figure  22  shows  a 
detailed  pressure  distribution  for  the  bow  region.  We  see  from  these  figures 
that  the  longitudinal  pressure  gradient  near  the  keel  is  favorable  and  later 
becomes  adverse.  The  pressure  gradient  in  transverse  direction  decreases 
rapidly  from  the  keel  to  a minimum  value  and  then  increases  continuously  up  to 
the  free  surface.  Due  to  this  rapid  pressure  variation  in  the  bow  region, 
preliminary  boundary-layer  calculations  showed  flow  separation  and  required  an 
approximate  procedure  to  generate  the  solutions  for  x < 22.5m.  After  that 
(x  > 22.5m),  the  three-dimensional  boundary-layer  calculations  were  performed 
for  a given  inviscid  pressure  distribution.  The  initial  conditions  at  x = 22. 
were  generated  by  solving  the  boundary-layer  equations  in  which  the  z-wise 


derivatives  for  a constant  z were 
neglected. 

Figures  23  to  25  show  some  of  the 

O 

computed  results  for  Rl  = 3 x 10°. 
Figure  23  shows  the  variation  of  cp, 
Cf,  Re,  Hu  and  ew  at  the  cross- 
planes of  x = 30m,  60m,  105m,  165m, 
210m.  Typical  streamwise  velocity 
profiles  at  x = 105m  and  z = 0.2 
are  shown  in  figure  24  and  typical 
crossflow  velocity  profiles  at  x = 

60m  are  shown  in  figure  25.  As  can 
be  seen  from  these  figures,  the  cross- 
flow  velocity  profiles  show  great 
variations  and  indicate  clearly  the 
flow  reversal  that  takes  place  in 

the  crossplanes.  This  implies  that 
differential  methods  based  on  two- 
dimensional  and/or  small  crossflow 


Figure  22.  Pressure  distribution  for 
the  bow  region  of  ship 
model  5350. 


i 
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approximations  as  well  as  methods 
based  on  Integral  methods  are  not 
adequate  to  boundary- layer  calcula- 
tions on  ship  hulls.  Other  interestlno 
results  that  emerge  from  these  calcu- 
lations are  the  sudden  Jumps  of  the 
limiting  crossflow  angle  from  positive 
to  neqatlve,  and  the  thickening  of  the 
boundary  layer  in  the  corner  region  of 
the  crossplanes.  The  jumps  of  the 
crossflow  angle  indicates  the  con- 
vergence of  the  flow  from  both  sides 
of  the  corner  region  and,  hence. 
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enhances  the  thickening  of  the  boundary 
layer.  This  thickening  of  the  boundary 
layer  in  the  corner  region  of  ship 
hulls  has  been  verified  experimentally 
by  Hoffman*’1. 


Figure  25.  Computed  crosswise  velocity 
profiles  for  ship  model 
5350,  for  Ri  * 3x10^ 
along  x * 60m  coordinate 
line. 


VI.  CONCLUDING  REMARKS  AND  FUTURE  WORK 


According  to  the  studies  presented  in  this  paper,  the  three-dimensional 

boundary  layers  on  ship  hulls  can  be  computed  very  efficiently  and  effectively. 

The  turbulent  model,  as  in  two-dimensional  flows,  again  yields  satisfactory 

results  for  three-dimensional  flows.  This  has  been  demonstrated  by  Soejima 
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and  Yamazaki  who  also  have  applied  the  present  turbulence  model  to  compute 
three-dimensional  boundary  layers  on  ship  hulls.  However,  there  are  additional 
studies  and  problem  areas  that  need  to  be  considered  and  investigated  before 
the  present  method  can  become  a more  effective  tool  to  design  ships.  They 
are  briefly  discussed  below. 

6.1  Generation  of  Initial  Conditions  on  Arbitrary  Bow  Configurations 

In  section  5-3,  we  presented  calculations  for  the  ship  model  5350  and 
mentioned  that  due  to  flow  separation  in  the  bow  region,  we  had  to  start  the 
boundary-layer  calculations  at  some  distance  away  from  the  bow.  Additional 
studies  are  required  to  generate  the  initial  conditions  on  the  bow.  These 
studies  can  lead  to  a better  design  of  bow  configurations  and  to  better 
handling  of  bilge  vortices,  which  contribute  to  the  total  drag  of  the  ship. 
However,  this  is  by  no  means  an  easy  task.  Consider,  for  example,  the  ship 
model  5350  discussed  earlier.  A sketch  of  the  bulbous  nose  with  a plausible 
inviscid  streamline  distribution  is  shown  in  figure  26.  We  may  assume  that 


Figure  26.  Pattern  of  streamlines  near  the  bow  of  ship  model  5350. 


the  ship  is  symmetrical  about  the  keel  plane  and  there  is  a nodal  attachment 
point  on  the  bulbous  nose  at  B.  If  the  ship  is  floating,  then  the  water 
line  is  determined  by  conditions  of  constant  pressure  and  zero  normal  velocity. 
Hence  the  intersection  of  the  plane  of  symmetry  with  the  water  line  at  A 
is  a saddle  point  with  the  streamlines  of  the  inviscid  flow  converging 
on  A along  the  line  BA  and  diverging  along  an  orthogonal  direction.  It 
is  known  that  the  boundary- layer  equations  can  always  be  solved  at  b but 
that  at  A the  situation  is  more  complicated  and  furthermore  it  is  still  not 
entirely  clear  what  their  role  is  in  relation  to  the  general  solution.  It  is 
likely,  however,  that  provided  no  reversed  flow  occurs  at  A in  the  component 
of  the  solution  along  the  direction  BA,  then  separation  can  be  avoided  along 
this  line  by  appropriate  choice  of  design.  Furthermore,  if  separation  does 
occur,  its  effect  may  be  limited.  The  recently  developed  Cebeci-Stewartson 
procedure  however,  can  be  applied  to  the  present  problem  but  there  are  some 
hurdles  to  be  overcome. 

Of  particular  difficulty  is  the  choice  of  coordinate  system  on  which  to 

compute  the  solution  and  to  join  it  with  the  already  well-established  method 

downstream  of  CD.  We  have  seen  that  in  the  case  of  the  prolate  spheroid 

20 

(see  Cebeci,  Khattab  and  Stewartson  ) it  is  helpful  to  have  a mesh  which  is 
effectively  Cartesian  near  the  nose  and  the  methods  which  were  used  to  produce 
it  in  the  earlier  study  are  applicable  to  any  body  which  can  be  represented  by 
a paraboloid  of  revolution  in  the  neighborhood  of  the  nose.  Now  here  we  have 
a paraboloid  near  B but  not  one  of  revolution,  but  we  believe  that  the 
necessary  generalization  is  possible.  The  mesh  now  has  to  match  with  that 
which  has  proved  convenient  downstream  of  CD.  Again  we  believe  that  a smooth 
transition  can  be  achieved  by  building  into  the  mesh  sides,  right  from  CBA,  an 
appropriate  spacing  such  that  the  points  of  a uniform  mesh  on  CD  are  also 
points  of  this  mesh  although  not,  of  course,  at  a constant  value  of  one  of  the 
coordinates.  Our  evidence  for  this  is  based  on  a successful  scheme  that  we 

have  already  worked  out  for  the  prolate  spheroid,  Cebeci,  Khattab  and 
20 

Stewartson  . 

Other  aspects  that  reed  further  study  include  the  condition  at  the  water- 
line section.  It  has  been  usual  to  assume  that  the  normal  velocity  is  zero  at 
the  undisturbed  free  surface.  This  is  not  quite  correct  and  the  error  may 
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have  implications  for  the  nature  of  the  solution  near  A and  especially  the 
question  of  separation  along  BA.  Even  if  separation  does  occur,  it  may  be 
possible  to  handle  the  post-separation  solution,  since  it  probably  extends 
only  over  a limited  region  of  the  ship,  by  means  of  an  interaction  theory, 
i.e.  modifying  the  inviscid  flow  by  means  of  a displacement  surface. 

6 . 2 Viscous-Inviscid  Flow  Interaction 

The  present  boundary-layer  calculations  are  done  for  a given  pressure 
distribution  obtained  from  an  inviscid  flow  theory.  In  regions  where  the 
boundary-layer  thickness  is  small,  the  inviscid  pressure  distribution  does 
not  differ  much  from  the  actual  one;  as  a result,  the  boundary-layer  calcula- 
tions are  satisfactory  and  agree  well  with  experiment,  see  for  example,  the 
papers  by  Cebeci,  Kaups  and  Moser^  and  by  Soejima  and  Yamazaki^9.  When  the 
boundary-layer  thickness  is  large,  which  is  the  case  near  the  stern  region, 
the  effect  of  viscous  flows  on  the  inviscid  pressure  distribution  must  be 
taken  into  account.  One  possible  way  this  can  be  done  is  to  compute  the  dis- 
placement surface  for  a given  inviscid  pressure  distribution  and  iterate. 

Such  a procedure  is  absolutely  necessary  to  account  for  the  thickening  of  the 

19 

boundary  layer  as  was  observed  by  Soejima  and  Yamazaki  . 

6.3  Prediction  of  Wake  Behind  Ship  Hulls 

The  present  boundary-layer  calculations  can  be  done  up  to  some  distance 
close  to  the  stern;  after  that,  flow  separation  occurs.  Since  one,  and 
probably  the  biggest,  reason  why  there  is  interest  in  boundary-layer  calcula- 
tions on  ship  hulls,  is  the  calculation  of  drag  of  the  hull,  additional 
studies  should  be  directed  to  perform  the  calculations  in  the  separated  region 
and  in  the  wake  behind  the  ship.  Recent  calculation  methods  developed  by 
Cebeci  and  his  associates  at  Douglas  for  two-dimensional  wake  flows  and  for 
separated  flows  by  using  inverse  boundary-layer  theory  are  appropriate  for 
these  purposes. 
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APPENDIX  A 

DESCRIPTION  OF  THE  BOUNDARY-LAYER  COMPUTER  PROGRAM 


This  appendix  describes  in  detail  the  computer  program  which  calculates 
three-dimensional  boundary  layers  on  arbitrary  ship  hulls  based  on  the  non- 
orthogonal  coordinate  system  and  the  numerical  method  described  in  the  text. 
The  computer  program  consists  of  two  separate  subprograms:  the  GEOMETRY  and 
BOUNDARY-LAYER  subprograms.  The  GEOMETRY  subprogram  requires  the  card  input 
of  the  ship  hull  geometry  section  by  section  and  the  corresponding  inviscid 
velocity  components  in  a Cartesian  coordinate  system.  It  performs  the 
following  functions:  (1)  smooths  the  input  data  by  the  use  of  Fourier  series 
expansion;  (2)  maps  each  section  defining  the  ship  hull  onto  a unit  circle 
and  defines  the  transverse  stations,  z,  for  boundary-layer  calculations; 

(3)  computes  the  geometric  parameters,  h-j , t^,  K] , l<2,  K-j 2 » K21 » S1 » and  9 
and  the  inviscid  velocity  components  along  the  boundary-layer  coordinate 
lines.  In  the  steps  2 and  3,  numerical  interpolation  and  differentiation  are 
frequently  encountered.  These  operations  are  accomplished  by  the  use  of  the 
cubic  spline  function.  Because  of  the  sensitivity  of  the  spline  function  to 
smoothness  of  data,  the  input  data  (which  are  sections  defining  the  ship  hull 
and  velocity  components)  have  to  be  fairly  smooth;  if  not, step  1,  which  is 
optional  in  the  program,  should  be  executed.  Since  the  mapping  of  an 
arbitrary  section  onto  a unit  circle  proceeds  clockwise,  it  is  convenient  to 
use  the  right-handed  Cartesian  coordinate  system  for  defining  the  ship  hull. 
Hence,  data  defining  each  section  must  be  input  clockwise  from  the  keel  to  the 
free  surface. 

The  boundary- layer  subprogram  requires  the  geometric  parameters,  hi, 
etc.,  and  the  inviscid  velocity  distributions  along  the  boundary-layer 
coordinate  lines,  which  are  the  output  of  the  geometry  subprogram  and  are 
stored  in  tape  unit  1.  In  addition,  physical  variables  and  program  control 
parameters  are  input  to  indicate  the  unit  Reynolds  number,  transition 
location,  and  specify  the  initial  condition  and  to  choose  the  grid  points 
across  the  boundary  layer.  Details  of  the  input  instruction  will  be  described 
later.  The  program  logic  and  structure  of  the  boundary- layer  subprogram  are 
rather  complicated;  the  basic  flow  chart  is  shown  in  Figure  A-l . 
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The  computer  program  is  written  In  Fortran  IV  for  the  IBM  370  system. 

When  other  systems  are  used  to  run  the  program,  It  Is  evident  that  slight 
modifications  may  be  required  to  meet  the  specific  requirements  of  that 
system.  Three  external  units  are  used  in  the  program  for  data  storage.  Unit  1 
(defined  as  TAPEDT)  Is  used  in  both  the  geometry  and  boundary-layer  programs 
for  storing  the  geometric  parameters  and  freestream  velocity  distributions 
generated  from  the  geometry  program.  Units  2 and  3 (defined  as  TAPEGP  and 
TAPEPF)  are  exclusively  used  in  the  boundary-layer  program  as  direct-access 
storage;  unit  2 stores  the  body  geometry  and  freestream  velocity  distribu- 
tions; and  unit  3 stores  the  pressure  parameters  and  boundary-layer  profiles. 
The  region  size  needed  to  execute  the  program  is  about  200K  bytes  based  on 
41  longitudinal  stations,  29  transverse  stations  and  61  grid  points  across 
the  boundary  layer. 

The  detailed  input  and  output  instructions  of  the  program  are  provided 
below. 

A-l  Input  for  the  Geometry  Program 

Card  1 contains  the  title  of  the  problem  under  consideration.  The  input 
is  punched  as  80-column  alphanumeric  field  as  shown  below. 


Load  Sheet  for  Card  1 

Card  2 requires  the  following  information  to  be  specified.  The  input  is 
punched  in  (415)  format. 


cc  1-5  NXTL  total  number  of  the  sections  defining  the  ship  hull, 

3 < NXTL  <41.  This  is  also  the  maximum  x-stations  (NX) 
in  boundary-layer  calculations. 
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cc  6-10  NZTL  maximum  transverse  stations  (NZ)  to  be  used  in  boundary- 

layer  calculations,  3 < NZTL  <_  29. 

cc  11-15  IM00TH  this  flag  controls  whether  the  input  data  (body  geometry 

and  freestream  velocity  distributions)  needs  to  be  smoothed 
= 0 no  smoothing  of  the  input  data 
= 1 smooth  the  input  data 


cc  16-20  I PRINT  this  flag  controls  printing  of  the  input  data 

= 0 no  printing  of  the  input  data 
= 1 print  the  input  data 


Card  3 specifies  the  following  information  related  to  the  ship  hull  sections. 
The  input  is  in  (F10.0,  315)  format. 
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Load  Sheet  for  Card  3 


cc  1-10  X(I) 
cc  11-15  INI 
cc  16-20  I QUAD 

cc  21-25  ICORN 


dimensional  x-value  of  Ith  cross  section  of  the  ship  hull, 
in  feet  or  meters. 

total  number  of  the  points  on  Ith  cross  section,  INI  <_  29. 
This  quantity  may  be  different  from  one  section  to  another. 

this  flag  identifies  whether  the  cross  section  is  completely 
submerged  under  the  water  surface  or  not 
= 1 the  cross  section  is  intercepted  by  the  water  surface 
= 2 the  cross  section  is  completely  submerged  under  the 
water  surface 

this  flag  indicates  whether  the  section  has  trail ing-edge 
corner  point  defined  in  Section  III 
= 0 there  is  no  corner  point  in  the  section 
= 1 there  are  corner  points  in  the  section 


Card  (Set)  4 contains  the  coordinates  of  the  section  and  corresponding  inviscid 
velocity  components  in  Cartesian  coordinate  systems.  The  input 
is  punched  in  (5F10.0)  format.  There  are  INI  cards  of  this  type 
in  each  section. 


1 1 ; 3 4 

• 

H | 

<■1  7 e 9|io 

1 Jl3 

14 

~n — i — 

15116  1711  20 

i 1 I 

J | t • 

2 5|?V*Ufc|2*rK 

j 1 [) ?l 3 .** 

j i i i i ' j r 

1 1 1:  - k 1 1. 

V 

B 

ZB 

UE 

VE 

WE 

i i i i i i i i_i_ 

-i  i i i i i i i i 

-iii  -i-i-  i i i-i. 

j 1 1 1.1. 

,,,,,,,,, 

Load  Sheet  for  Card  4 


cc  1-10  YB 


dimensional  y-value  of  the  cross  section,  in  feet  or  meters. 


cc 

11-20 

ZB 

cc 

21-30 

UE 

cc 

31-40 

VE 

cc 

41-50 

WE 

dimensional  z-value  of  the  cross  section,  in  feet  or  meters. 

nondimensional  inviscid  velocity  component  in  x-direction, 
ue/Uref. 

nondimensional  inviscid  velocity  component  in  y-direction, 
Ve^^ref • 

nondimensional  inviscid  velocity  component  in  z-direction, 

we/uref • 


It  should  be  noted  that  Card  3 and  Card  (Set)  4 have  to  be  input  for  all 
cross  sections  (NXTL) 


A-2  Input  for  the  Boundary-Layer  Program 

Card  1 contains  the  title  of  the  problem  under  consideration.  The  input  is 
punched  as  80-column  alphanumeric  field  as  shown  below. 
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Load  Sheet  for  Card  1 


Card  2 contains  the  following  information  to  be  specified.  It  is  punched 
in  (615)  format. 
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Load  Sheet  for  Card  2 


cc 

1-5 

NXT 

cc 

6-10 

NZT 

cc 

11-15 

NXSTRT 

cc 

16-20 

NTR 

the  number  of  the  last  x-station  (NX)  to  be  calculated, 

3 < NXT  < NXTL 

the  number  of  the  last  z-station  (NZ)  to  be  calculated, 

3 < NZT  < NZTL 

the  number  of  x-station  (NX),  referred  to  the  geometry 
data,  where  boundary-layer  calculations  begin. 

the  number  of  x-station  (NX)  where  turbulent  flow  calcu- 
lations begin.  For  laminar  flow  only,  set  NTR  > NXT; 
for  turbulent  flow  only,  set  NTR  = 1.  NTR  will  be  over- 
ridden if  laminar  separation  occurs. 
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cc  21-25  NSWCH 


the  number  of  x-station  (NX)  counted  from  NXSTRT  where 
switching  to  "zigzag"  scheme  is  made.  If  no  switch  is 
desired,  set  NSWCH  > NXT.  If  NTR  = 1,  NSWCH  may  be 
arbitrary  depending  on  the  transverse  pressure  gradient; 
if  NTR  f 1,  NSWCH  ^3  is  recommended. 

cc  26-30  JPRINT  point  interval  by  which  the  profiles  are  printed  out.  If 

JPRINT  = 1,  the  velocity  profiles  are  completely  printed 
out  from  the  wall  to  the  boundary-layer  edge. 

Card  3 contains  the  following  information  to  be  specified.  The  input  is  in 
(5F10.0)  format. 


■MMBffl 

BBBBBBBBBB 

B3SBEB99S9 

BBMBBBBES 

ETAE 

VGP 

DETA(l) 

CNU 

UREF 

i i i i i i i i i 

1 1 1 1-1  1 1 I -L- 

i i i i i i-  i i -i  j 

it  l 1 1 1 l- 

i i i i i — i — i-.-j._i— 

Load  Sheet  for  Card  3 


cc  1-10  ETAE 


cc  11-20  VGP 


cc  21-30  DETA(l) 


transformed  boundary-layer  thickness,  for  the  first 
station.  A value  of  8 is  usually  sufficient.  For  NX  > 1 , 
it  is  computed  internally. 

variable  grid  parameter,  K >_  1.0.  For  laminar  flow, 

K = 1.0  is  sufficient.  For  turbulent  flow,  it  is  a 
function  of  the  Reynolds  number,  R 


l-i  spacing  at  the  wall.  For  laminar  flow  only, 
is  suggested.  For  turbulent  flow,  it  is  also 


initial  Am 
Ani  = 0.2 

a function  of  R^.  The  suggested  K and  Am  values 
for  different  Reynolds  number  are  provided  below. 


R 

00 

106 

107 

108 

109 

K 

1.10 

1.15 

1.20 

1.25 

An-| 

0.015 

0.01 

0.0075 

0.0050 

CC  31-40  CNU 

kinematic 

viscosity 

of  the 

fluid,  v 

, in  ft 

cc  41-50  UREF  reference  velocity,  uref,  in  ft/sec  or  m/sec. 


A-3  Output  for  the  Geometry  Program 

The  output  of  the  geometry  program  includes  the  printout  of  the  input 
data,  the  computed  geometrical  parameters,  and  the  inviscid  velocity  components 
along  the  boundary- layer  coordinate  lines.  The  latter  are  also  stored  in 
external  store  unit  (called  TAPEDT)  and  serve  as  a part  of  the  input  data  to 
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Li 


; 


the  boundary- layer  program.  The  notation  used  in  printing  out  the  input  data 
is  the  same  as  that  used  in  the  input.  Therefore,  only  the  meaning  of  the 
calculated  geometrical  parameters  is  Interpreted  here. 

NX  station  number  of  the  x-station 

x x-value  at  NX  station 

z normalized  transverse  coordinates  in  nonorthogonal  boundary-layer 

coordinate  system,  measured  from  the  keel  (z  = 0.1)  to  the  free 
surface  (z  = 1.0) 

YB  input  or  interpreted  y-value  on  the  ship  hull 

ZB  input  or  interpreted  z-value  on  the  ship  hull 

HI  metric  coefficient  hi  associated  with  the  coordinate  x 

H2  metric  coefficient  associated  with  the  coordinate  z 

K1  geodesic  curvature,  K-j , of  the  x-coordinate  line 

K2  geodesic  curvature,  l<2,  of  the  z-coordinate  line 

SI  the  physical  surface  distance  Si  measured  along  the  x-coordinate 

line  from  the  first  x-station 

K12  geometric  parameter, 

K21  geometric  parameter,  1^ 

THETA  angle,  e,  between  the  coordinate  lines  in  radian 

UE  longitudinal  velocity  component,  ue,  in  the  boundary-layer 

coordinates 

WE  transverse  velocity  component,  wg,  in  the  boundary-layer 

coordinates 

2 

Cp  pressure  coefficient,  2(p  - Pm)/puoo 

GIRTH  normalized  surface  distance  along  z-coordinate  line,  measured 
from  the  keel  to  the  free  surface. 

A-4  Output  for  the  Boundary-Layer  Program 

The  output  of  the  program  includes  printout  of  the  geometrical  data 
passed  from  the  geometry  program  on  external  storage  unit  1,  as  well  as 
tables  of  boundary-layer  profiles  and  some  important  boundary-layer  parameters. 
The  notation  for  the  geometrical  data  is  the  same  as  that  in  the  geometry 


program;  therefore,  only  those  associated  with  the  boundary-layer  parameters 
are  given  below. 


1 . The  Profile  Data 

J point  number:  profiles  are  printed  from  the  wall  outward.  Not 
all  points  are  printed  and  the  interval  is  controlled  through 
the  input  variable,  JPRINT. 

ETA  nondimensional  boundary-layer  variable  n. 

Y normal  distance  from  the  wall,  y,  ft  or  m. 

U nondimensional  velocity  in  x-direction  in  the  boundary  layer, 
f*  = u/ue. 

W nondimensional  velocity  in  z-direction  in  the  boundary  layer, 
g*  = w/uref. 

V derivative  of  U with  respect  to  n,  f". 

T derivative  of  W with  respect  to  n,  g". 

US  velocity  component  in  the  boundary  layer  along  the  inviscid 
streamline  direction,  u$. 

UN  velocity  component  in  the  boundary  layer  normal  to  the  inviscid 
streamline  direction,  un. 

USE  total  velocity  at  the  edge  of  the  boundary  layer,  u$e. 

BETA  crossflow  angle  in  degree,  8,  defined  as 

8 = tan_1un/us 

B nondimensional  total  eddy  viscosity,  b = 1 + e/v 

2.  Boundary-Layer  Parameters 

CFS  local  skin-friction  coefficient  in  the  inviscid  streamline 
direction,  cf  = 2tu  /puf 

a S e 

CFN  local  skin-friction  coefficient  normal  to  the  inviscid 
streamline  direction,  Cfn  = 2tw  /pufg 

DLSTS  displacement  thickness,  6*,  in  feet  or  meters,  defined  as 

oo 

5* -f  ^ “ us/use)dy 


i i 
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lllTAS  momentum  thickness,  o^,  In  feet  or  meters,  defined  .is 

°11  " I us/%  0 “ us/useRv 

MS  shape  factor, 

RTHTA  Reynolds  number  based  on  momentum  thickness,  R(1  - i']|Us  /v 
Cf X local  skfn-frlctlon  coefficient  In  x-dlrectlon. 


fx 


-/ pU*  f 

x ref 


KL Y Reynolds  number,  Rx  * useS./v 
Cf'XB  average  skin-friction  In  x-dlrectlon, 

cf-  * /V-  dA/A 

fx  J tx 

ARIA  total  ship  hull  surface  area  measured  from  the  starting  station 
to  the  current,  station 


mi  ia 


boundary-layer  thickness,  i1'.  In  feet  or  meters 
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